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ABSTRACT 



Tidal disruption of stars by massive black fioles produce transient accretion 
flows that flare at optical, UV, and X-ray wavelengths. At late times, these 
accretion flows may launch relativistic jets that can be detected through the 
interaction of the jet with the dense interstellar medium of the galaxy. We present 
an upper limit for the flux density of a radio counterpart to a tidal disruption 
event detected by GALEX that is a factor of 6 below theoretical predictions. We 
also examine existing radio surveys for transients with a time scale of 1 year and 
use these to set a 2a upper hmit on the rate of tidal disruption events producing 
relativistic jets of ~ 14 x 10^^ Mpc^'^ y^^. This rate is an order of magnitude 
lower than the highest values from theoretical models and is consistent with 
detection rates from optical and X-ray surveys. 



Subject headings: radio continuum: general 
uum: galaxies — surveys 



galaxies: active — radio contin- 



1. Introduction 



Massive black holes at the centers of galaxies disrupt stars that pa ss within t he tidal ra- 



dius of the black hole, resulting in a process of transient accretion (e.g. . iRees 119881). Dynami- 
cal m odels predict that these processes occur every 10^ to 10^ y per galaxy (IWang fc Merritt 



2004 ). The accretion process leads to a flare that peaks at X-ray or UV wavelengths and has 
a power-law decline with a characteristic timescale of months. Such events are recognized 
by the flare and th e transition from a normal galaxy to a state with an AGN spectrum. 
Gezari et al.l (|2009[ ) summarize the current state of our knowledge regarding tidal disruption 
events detected at X-ray, UV, and optical wavelengths. In total, 10 events have been char- 
acterized. Synoptic optical surveys may detect tens of such events and prove powerful in 
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characterizing tidal disruptio n events and determining their statistics ( IStrubbe fc Quataert 
2OI0I : Ivan Velzen et al.ll201ol ). 



A significant astrophysical problem is the coupling between accretion and the launching 
of jets. Such systems are ubiquitous yet there is no complete theory that can predict jet 
properties based on a known accretion flow. Tidal disruption events represent a unique 
laboratory for exploring the accretion-outflow coupling under a range of different accretion 
conditions over a short period of time. If accretion leads to a relativistic outflow, this 
will likely radiate synchrotron emission, which can be detected at radio wavelengths. A 
recent model explores r adio emission that results from the interaction of the outflow with 
the interstellar medium ( IGiannios &: Metzgerll201ll ). Radio surveys have the advantage that 
they are unaffected by extinction in the dense nuclear environments of the host galaxy; 
however, there are no known detections of radio emission from tidal disruption events. 

In this paper we examine what constraints that existing and future radio surveys can 
place on tidal disruption events. In § [21 we briefly summarize a new model for radio emission 
from tidal disruptions. In § [31 we discuss the relevant survey data and identify candidate 
events. In § [H we use these results to set constraints on the rate of tidal disruptions that 
produce relativistic jets. In § [5l we present archival Very Large Array results for a tidal 
disruption event detected by GALEX. In § [6l we summarize our results. 



2. Radio Emission from Jet-Induced Reverse Shocks 



Giannios fc Metzgerl (120 111 , hereafter, GM) introduce a model in which a reverse shock 
from the interaction of a jet with the dense interstellar medium of the galaxy produces an 
afterglow. The afterglow has a peak flux density 

where D is the distance to the source and C is a constant of order unity that depends on 
details of the jet and ISM physics. The spectrum peaks at frequency of ~ 25 GHz and 
extends as v^^^ to lower frequencies before reaching the self-absorption frequency. As the 
source cools, the peak flux density propagates to lower frequencies with proportionality z/^/^. 
The reverse shock spectrum is expected to peak ~ 1 year following the tidal disruption event 
and have a duration of 1 year or less. This afterglow has the important property that it emits 
isotropically. Accordingly, events are not limited by the beaming angle of the relativistic jet 
that creates the relativistic shock. 



A significant challenge in identification of a tidal disruption event in a nearby galaxy 
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may be confusion with galactic emission. In this case, the tidal disruption would not appear 
as a true transient but rather as significant variability in a detected source. The amplitude 
of this variability must be large in the model of GM, however. For example, the galaxy NGC 
891 has a fiux density of 700 mJy at 1.4 GHz. If moved from its current distance of 7 Mpc 
to 1 Gpc, it would have a fiux density of 35 microJy, a factor of 10 below the predictions of 
GM for a tidal fiare at this distance. Thus, in practice tidal disruption events will appear as 
transients or extreme variables in most surveys. 



3. Survey Results 

We summarize here existing radio surveys that have sensitivity on timescales of ~ 1 year 
and are, therefore, sensitive to tidal disruption events. In Table [H we give the frequency 
(/), limiting fiux density (Fum), survey area number of epochs (iVg), and number of 
candidate events (A^can)- The number of epochs is the number of independent year- long 
observations. 



The archival VLA survey (iBower et al.l 120071 ) at 5.0 and 8.4 GH z observed a single 



pointing weekly for 20 years, producing nearly 1000 independent epochs. iBower et al.l (120071 ) 
produced maps and searched for transients on time scales of a single epoch, 2 month integra- 
tions, and 1 year integrations. No radio transients were detected from the 1-year integrations. 
Two radio transients were detected from the 2-month integrations but neither is associated 
with the nucleus of an optically-detected galaxy. A characteristic sensitivity of 0.1 mJy is 
chosen to represent the overall survey at 2-month and 1-year sensitivity. At this fiux density, 
the effective survey area over all epochs is 0.26 deg^ at 1 year. This quantity is equal to 
N^fl. Dividing over A'^e = 20, we find the per epoch area of 0.013 deg^. Limits for two-month 
integrations are similar to the 1-year results. Angular resolution of the survey ranged from 
< l"to ~ 15". 

An analys i s of 3 C 286 observations from the VLA archives at 1.4 GHz was presented by 



Bower &: Saull (1201 ll ). This paper explored daily variations in data obtained over 23 years 
in 1852 epochs. Since integrations longer than a single day were not produced, the single 
epoch fiux density threshold sets the transient sensitivity limit. This limit is relatively poor, 
limiting the overall sensitivity of this result to month- or year-long events. These data were 
taken in the C and D configurations of the VLA, giving characteristic angular resolution 
> 15". No radio transients were found above a characteristic fiux density of 70 mJy. The 
solid angle is set by the field of view imaged. 



The PiGSS-I (IBower et al, 



20101) and ATATS-I (IBower et al 



2010|) results are from Allen 



-4- 



Telescope Array (IWelch et al.l 120091 ) single-epoch surveys at 3.1 and 1.4 GHz, res pectively. 
These survey results are compared against the NVSS catalog (ICondon et al.lll998[ ) to deter- 
mine the rate of transient sources. The characteristic timescale of both of these efforts is 
any timescale longer than the integration timescale of NVSS, which is less than 1 minute. 
NVSS observations were conducted more than 1 decade earlier than the ATA observations. 
The PiGSS-I catalog is constructed from observations spread over 4 months in 2009 and has 
a resolution of 100". The ATATS-I catalog is constructed from observations spread over 4 
months in 2008 and has a resolution of ~ 200". 

Light curves of sources detected in the VLA, 3C 286, PiGSS-I, and ATATS-I surveys 
indicate no variability greater than a factor of a few in all surveys on long timescales. 



Levinson et a. 



(I2OO2I) and Gal- Yam et al.l ( 120061 ) presented analyses of a comparison of 



the VLA FIRST f lBecker et al.lll995f ) and NVSS 1.4-GHz catalogs in search of radio tran- 
sients. The characteristic separation time between these catalogs is ~ 1 year. A radio 
supernova was the only transient identified. The characteristic resolution of the two surveys 
was 5 and 45", respectively. The authors report no evidence of extreme variability in sources 
in either epoch. 



Bannister et al.l (|2010[ ) searched the archives of the Molongolo Synthesis Telescope to 



identify all fields observed multiple times over a 22-year period. The total surface area sur- 
veyed was 2776 deg^ at a flux density threshold of 6 mJy. All MOST observations were 
conducted at 0.843 GHz with a resolution of 45". The heterogeneous nature of the sample 
makes determination of the sensitivity of the survey at year and longer timescales difficult. 
Examination of Figure 5 in their paper indicates that approximately two-thirds of the ob- 
servations were paired with others on timescales greater than 1 year. Accordingly, we scale 
the total area by a factor of two-thirds. 

The MOST survey identified 53 variable and 15 transient sources. Of the 15 transient 
sources from MOST, only three are unambiguously identified with known objects (SN 1987A, 
Nova Muscae 1991, and GRO 1655-40). Another source is offset from the nucleus of a nearby 
spiral galaxy and, therefore, unlikely to be a tidal disruption event. We can also exclude 
three sources showing variability on timescales much less than 1 year. We further exclude 
four sources that exhibit multiple detections with intervening non-detections. Two of the 
remaining four sources have optical detections in archival data indicative of a quasar. The 
dates of these optical observations are unclear but the surveys from which the data are drawn 
began and ended several years prior to the transient detection epoch for both sources. The 
authors argue that these are scintillating AGN. The remaining two sources (J064149-371706 
and SUMSS J102641-333615) have a single detection with MOST, a non-detection with 
MOST about a year after the original detection, and no identified optical counterpart. We 
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note that both sources have apparent NVSS counterpart that the authors did not identify, 
indicating that these are variable sources rather than true transients. Thus, none of the 
MOST transients quahfy as plausible tidal disruption candidates. 

Additionally, two variable sources (J201524-395949 and J200936-554236) from MOST 
are possibly associated with the nuclei of spiral galaxies at redshifts of 0.02 and 0.032, 
respectively. J201524-395949 exhibits an increasing flux density over 10 years with amplitude 
of ~ 20 mJy and a factor of 4 total increase. J200936-554236 exhibits a 5 mJy rise in flux 
density over three years (corresponding to less than factor of 2 increase in flux density). 
The amplitude of variability in both sources falls well below the predictions of GM (~ 200 
mJy at this distance). Nevertheless, these sources broadly follow the expectations for tidal 
disruptions and are plausible candidates. 



4. Limit on the Tidal Disruption Rate Producing Jets 

We also tabulate calculated values for the limiting distance (Aim) at which tidal disrup- 
tions would be visible under the model of GM and an upper limit to the tidal disruption rate 
producing jets {ttd)- is computed directly from the limiting sensitivity of the surveys 
and explicitly making use of the decreasing peak flux density at lower frequencies (oc z/^/^). 
We compute the rate from 

-"can /r,N 

TTO = (2) 

where V is the volume surveyed. For the case of no candidates, we set the upper limit at 2a 
with Ncan = 3. We give ttd in units of 10~^ y~^ Mpc~^. 

These estimates clearly indicate that three surveys set the strongest limits. The VLA 
archival survey sets limits of r^D < 29 for both two-month and 1-year integrations, respec- 
tively. The FIRST-NVSS comparison sets an upper limit ttd < 14 while the MOST survey 
gives a rate limit vtd ^ 17. Note that in the case of MOST, the existence of plausible 
but not convincing candidates still leaves this value as an upper limit. Summing over all 
independent surveys, we estimate ttd ~ 4.3. We can place an upper limit of ttd < 14 at 
2a confidence. 



The radio upper limit is consistent with estima tes based on ROSAT (jPonley et al. 



2OO2I ) and SDSS observations (Ivan Velzen et al.l I2OIOI ). The latter estimates a rate N 



3 X 10^^ y~^ galaxy"^. For a space density of galaxies with massive black holes of ~ 3 x 
10"'^ Mpc~'^, the optical rate implies ttd ~ 1-0. The limit that we have set falls more than an 
order of magnitude below the most opti mistic model for the ev ent rate of tidal disruptions. 



which estimates rrn ~ 10 ^ y ^ Mpc ^ (jWang fc Merrittll2004j ). 
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Tidal disruption rates for radio detected events are likely to be lower than in the optical 
given that only a fraction of accreting systems typically produce radio jets. Thus, the radio 
limit can be considered consistent with the most optimistic optical constraints if < 10% of 
the systems produce relativistic jets. Similarly, changes in the jet physics model of GM can 
strongly influence the event rate limit. For instance, the rate ttd scales as the jet luminosity 
to the -3/2 power. 

Future surveys will be capable of significantly improving sensitivity. As GM note, the 
complete PiGSS results have the capability to detect tens of events. Non-detections would 
place significant constraints on ttd and the reverse shock model. Similarly, further analysis 
of data from the VLA archives can provide significantly more sensitivity. Year- or month- 
long integrations of the 3C 286 data would likely set sensitivity levels at 1 mJy or better 
while covering comparable area. The resulting sensitivity to tj-d would be in the range 
of 10. Use of additional sources or data at other frequencies could lead to an order of 
magnitude improvement. The VAST intermediate depth survey with ASKAP will achieve 
250 micro Jansky sensitivity over 10^ deg^ in multip le passes and have substantial sensitivity 
to phenomena of this kind (jChatterjee et al.ll2010l ). 



5. VLA Images in the Field of GALEX J141929+525206 



Gezari et al.l ( 120061 ) reported the detection of a UV flare from an early-type galaxy with 
no prior evidence for AGN activity at z=0.37. The source was undetected in 2003 June 21 
- 29, first appeared in 2004 March 25 - June 24, and then decayed steadily until it was not 
detected in 2006 March 5 - 7. 

We used archival VLA data of observations in the field of the tidal flare candidate source. 
Data were obtained at a frequency of 1.4 GHz in spectral line mode to enable full field of 
view imaging on 2003 December 2, 2 005 April 18, and 2005 June 20 as part of a campaign 



to survey the extended Groth strip (llvison et al.ll2007| ). No source was detected in these 



individual epochs at the location of the GALEX source with image rms of < 50 ^Jy. We 
integrate together all data to produce an image with 15 /uJy rms and a synthesized beam 
4.0" X 3.9" (Figured]). Again, no source is detected in the image of the GALEX source. A 
bright, extended source in the field is associated with the galaxy SDSS J141930. 12+525159. 

The model of GM predicts a radio flux of 0.2 mJy for this source at a frequency of 1.4 
GHz. This is a factor of 6 higher than our 2a upper limits in the integrated image. Thus, 
our observations, which are coincident with the first appearance of the UV flare and extend 
to 1.5 y after the peak, demonstrate that this source likely did not produce a relativistic jet 
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that led to a reverse shock in the interstellar medium, or that the reverse shock emission was 
relatively short-lived. 



In addition, we note that Ivan Velzen et al.l (l2010f ) obtained upper limits for a radio 
counterpart of ~ 0.1 mJy at 8.5 GHz to tidal disruption event TDE2 in a galaxy at a 
distance of 1.05 Gpc. The flux densities fall well below the GM peak-flux predictions of 
~ 1.1 mJy. These limits were obtained 7 and 92 days following the first initial optical 
detection, and so may not probe the late-time shock event. 



6. Summary 

We have presented here a summary of existing radio surveys with sensitivity to transients 
with timescales of ~ 1 y. These surveys provide the first constraint based on radio data of the 
rate of tidal disruption events around massive black holes in the nearby Universe. This upper 
limit is consistent with rates set via optical methods. Note that the rate we have set depends 
on details of the GM model for the jet and its interaction with the interstellar medium of the 
galaxy. Factors such as the fraction of systems that produce jets, the efficiency with which 
accretion power is converted to jet luminosity, and the density of the interstellar medium 
can significantly affect the estimated rates. The non-detection of a radio counterpart to a 
tidal disruption event demonstrates the degree to which we do not understand the accretion- 
outflow characteristics of systems of this kind. 

Ongoing surveys such as PiGSS, future surveys such as VAST, and analysis of archival 
data can contribute to our constraints on radio properties of tidal disruptions. Equally 
important is follow-up of candidate events with sensitive observations using the EVLA. 
Unambiguous detection of a single radio counterpart would provide significant insight. Future 
optical surveys such as those carr ied out by Pan-STARRS are predicted to detect 0(10) 



events per year (iGezari et al.ll2009l ). Following up all events at a sensitivity 10 times deeper 
than the predictions of GM would provide a window into what mechanism in the accretion- 
outfiow coupling leads to the paucity of radio emission from tidal disruption events. 
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Fig. 1.— Image of the field of the tidal disruption event GALEX J141929+525206. The 
cross marks the position of the GALEX source. Contours are -3, 3, 4.2, 6, 8.4, 12, 16.8, 
and 24 times the rms noise of 15 //Jy. The grey scale ranges from -30 /^Jy to 150 /^Jy. The 
synthesized beam is shown in the lower left. 
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Table 1. Surveys with Long Timescale Sensitivity 



Name 


/ 








^ ^ can 




rTD 




(GHz) 


(mJy) 


(sr) 






(Gpc) 


(10-^ y-^ Mpc-3) 


VLA 


5.0 


0.1 


3.9C-06 


20 





2.99 


< 29 


3C286 


1.4 


70.0 


1.6e-04 


23 





0.08 


< 29042 


PiGSS-I 


3.1 


2.0 


3.0e-03 


1 





0.59 


< 96 


ATATS-I 


1.4 


230.0 


2.1e-01 


1 





0.05 


< 3113 


MOST 


0.8 


14.0 


5.5^01 


1 


2 


0.16 


< 17 


FIRST-NVSS 


1.4 


6.0 


1.9e-01 


1 





0.28 


< 14 



